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Abstract 

We study a minimally coupled charged scalar field in a charged Lifshitz background. For z — 2, we find 
an analytic expression for the corresponding low energy retarded Green's function. Unlike the RN-AdS 
case, the position of the superfluid surfaces depends on the charge of the scalar field only through the IR 
scaling dimension. We show that by increasing the dynamical exponent, the dual theory becomes more 
stable. We also show that the background could suffer from an instability of the IR geometry leading to 
a bifurcating critical point. It also allows the existence of scalar hair, causing hybridized critical point. 
We have investigated stable and unstable regions in the parameter space. 



1 Introduction 



The celebrated gauge/gravity duality [IH3] provides a useful tool to study strongly coupled quantum field 
theories via classical gravity. Over the past few years, this approach has been applied to study certain areas 
in condensed matter physics. In this context, one of the most challenging problems is to find a framework to 
study strongly correlated systems, such as non- Fermi liquids (for excellent reviews see [H[S]). An interesting 
feature of such systems is their behavior near quantum critical points (QCPs). 

At the critical point, a physical system typically exhibits scaling symmetries which in general may be 
given by 

t — > \ z t, Xi^-Xxi, (1.1) 

where z is the dynamical exponent. Such a symmetry, known as Lifshitz scaling symmetry, is widely studied 
in condensed matter physics starting with [6], The case of z = 1 reduces to the relativistic case, while z > 1 
refers to anisotropic scaling symmetries. 

To holographically study theories with this scaling symmetry, one may consider a classical gravity on the 
following background [7^1 

, 2 r i( dt 2 dx 2 dr 2 \ , „, 

ds = L {-—*+ — + — (1-2) 

where L is the radius of curvature. This metric is invariant under (jl.ip together with r — > A r. Also 
analytic solutions describing black-holes and black-branes in asymptotically Lifshitz spacetimes is constructed 
by considering Einstein-Maxwell-dilaton theories [10] . 

Following our previous work [11] . where asymptotic Lifshitz solutions were probed by fermions, in this 
paper we will probe a charged Lifshitz black-brane by a charged scalar field in order to study the dual theory 
at the low energy limit. In this context, certain features of critical behavior, including bifurcating, hybridized, 
and marginal QCPs have been studied on RN-AdS background in four and five dimensions [12]|T3]. Unlike 
thermal phase transitions, these quantum phase transitions could not be fully described by the Landau- 
Ginsburg paradigm. We will investigate these behaviors in a non-relativistic model. 

We show that although this non-relativistic theory flows to an IR fixed point, the non-relativistic nature 
affects the critical behavior of the system via the dynamical exponent. Increasing the dynamical exponent 
suppresses the instabilities due to the violation of the BF bound in the near horizon AdS2 region. We will 
focus on the case of z = 2, which is claimed to be the margin of different physical behavior 14J. In this paper 
we will solve a scalar field analytically and work out the bosonic retarded Green's function at low frequency 
limit for standard quantization. The retarded Green's function has several zero modes, which could be a 
source of instability in the dual theory. The zero modes are located at 

1 + A- 
2 Vv k = -n, 

where n is a non-negative integer, is the IR scaling dimension, and A_ is the scaling dimension of the 
dual operator defined by 

v k = - \/'Sm 2 + 3fc 2 + 9 - 4e 2 , A_ = -2 + \/4 + m 2 . (1.3) 
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Unlike the case of RN-AdSs, the locus of zero modes does not depend on the charge of the scalar field 



As it has been mentioned in [7], although the metric is non-singular, it is not geodesically complete and has divergent tidal 
forces in the far infra-red. The situation where this divergence can be resolved is studied in [8]|9]. 
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explicitly [P5] . 

We will proceed by studying different types of phase transitions in this model. The first type of phase 
transition occurs because of the violation of the BF bound deep in the space-time geometry. This point is 
called the "bifurcating" critical point. At such a critical point, the static susceptibility is not divergent but 
bifurcates into the complex plane. We will show that our model supports such a behavior and investigate 
different aspects of the phase diagram. 

The second type of phase transition occurs when the system supports zero modes and leads to a hairy 
black-brane geometry. The onset of this behavior is known as the "hybridized" critical point. This critical 
point could be approached by turning on a double trace deformation. We study different regimes in the 
parameter space which support the hybridized critical point. We also show that the model could be de- 
formed by a relevant or an irrelevant double trace operator, corresponding to the alternative and standard 
quantizations. 

Approaching bifurcating and hybridized critical points at the same time leads to a third type of critical 
behavior. At such a point known as the "marginal" quantum critical point, we will study the Green's function 
which shows the behavior of marginal Fermi liquids. 

This paper is organized as follows. In Sec. 2 we briefly review the charged Lifshitz black-brane solutions 
and their near horizon geometries. In Sec. 3 we shall discuss the non-relativistic behavior of the IR scaling 
dimension. Also the low energy retarded Green's function of a charged scalar field is calculated analytically 
for the case of z = 2. In Sec. 4 we will study different types of critical behavior of the model. The last 
section is devoted to conclusions. The appendix contains some details of calculations. 



2 Charged Lifshitz black-brane 

In order to study strongly correlated systems at finite density, one may consider a charged Lifshitz black-brane 
as a possible gravity dual. In the case where the temperature is much smaller than the chemical potential, 
the extremal limit of the black-brane must be considered. Thus we will consider an extremal charged Lifshitz 
black-brane in the bulk. 

Pure Einstein gravity does not support Lifshitz geometry as a solution. One can obtain Lifshitz geometry 
either by coupling a massive gauge field, or by adding higher derivative terms to Einstein gravity [TB] . Asymp- 
totic Lifshitz black-branes could be obtained as solutions of Einstein-Maxwell-dilaton theorieo Consider the 
following action [TO] 
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16ttG 4 



R-2A- i(90) 2 - i e A ^(F«) 2 - l e x ^(F^) 2 



(2.1) 



This model admits a charged black-brane solution where z plays the role of critical exponent and A, Ai and 
A2 are given by 

(z + l)(z + 2) 



A = -- 7^2 "> A i = ~^^' A 2 = % /^1. (2.2) 

ZL- v z — 1 



2 Suporfluids and superconductors in which the gravitational theory includes a dilatonic field and are asymptotically AdS are 
studied in |17j . 
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with the following solutioij§ 

dr 2 1+r 2(z+1) r 2(z+1) 

* 2 = -r 2 *fdt> + %- f + / = 1 - 1±^_ + ^ 

4zrf +1 > 

4 1 )=- M ( 1) (l-r 2 ^), Af) = ^(l-^), (2.3) 
where the mass and charge of the black-brane is determined in terms of r$ and k, and 



2(z+l) 



^-M*-V[ ^ „00 = *SL_ (2 4) 



' o 



bmce diverges at the boundary, it can not be treated as a chemical potential of boundary theory and 
we are forced to couple matter fields to A^ 2 \ 
The Hawking temperature reads 



z + 2 



Using the general idea of AdS/CFT correspondence, the physics at low energy is governed by the near horizon 
modes. Thus we will study a charged scalar on the near horizon background. At zero temperature where 
r 2( z +i) = £±2 and j = ! _ 2(1 + z -iy-{z+2) + (]_ + 2z- 1 )r-^ z+1 \ setting 

r-l = ? * , - u , t=-r, (2.6) 
(z + l)(z + 2)£ e 

the near horizon background can be obtained by taking the limit e — > where one finds 



2 



ds 2 = ^ (-dr 2 + d£ 2 ) + *c 2 (2-7) 



and l z 2 — (z + l)(z + 2). In this limit the dilaton and gauge fields reads 



e v 



(i) 



4(1) - £ 4(2) _ 

4(z + 2)' T + T + 1 j 



The near horizon geometry takes the form of AdS2 x R and is invariant under the following scale transfor- 
mations 

t — y At, £ — ;> A£, Xj -> Xj. (2.9) 

Therefore the low energy physics is described by an emergent IR CFT. As a result, we would expect the 
model to exhibit superfluid surfaces whose physics is governed by an IR fixed point. The R 2 factor of the 
near horizon geometry together with the AdS2 factor leads to the so-called semi-local quantum liquid (SLQL) 
phase [18] . 



■i . a) a) 

Note that we have shifted the gauge fields by constants to make sure that g^AyAl ' remains finite. We have set L = 1 
and also with a proper choice of the parameters the radius of horizon has also been set to unity. 
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Figure 1: Left plot: Square of v k as function of dynamical exponent is shown for different values of charge. 
As the dynamical exponent increases, the systems tends to IR stable region. The absolute value of charge 
compensates reaching the stable region. Right plot: Square of v k as function of charge for different values 
of z. As the charge increases, v k tends to zero and finally takes pure imaginary values which causes IR 
instability. Notice that for larger dynamical exponents the systems reaches the unstable region softer. In 
both plots to = 1 and k = 0. 



3 Analytic Green's function 



In order to find the retarded Green's function one must solve the scalar field equation of motion on background 
(|2.3I) . which reads 



r' z - x d T (r z+a fd r cl>) + 



(w + qA 



(2)N2 



r 2z f 



= 0. 



(3.1) 



In the low frequency limit, since / has a double zero at the horizon, one can not take the w — >■ limit strait 
forwardly. To treat this regime, following |18j . the geometry may be divided into IR and UV regions. In the 
IR region one can not neglect the u> term due to the singularity of / near the horizon. But in the UV region 
where / is regular, one can set to — at first order. Matching the solutions in the intersection region could 
lead to the retarded Green's function of the whole geometry at leading order in ui. 



3.1 Low energy behavior and IR Green's function 

Solving Eq. (|3.1[) in the IR geometry (|2.7p leads to the following IR Green's function (for details see the 
appendix) 

r t,,\ a -inu k r(-2i/ fc ) r(± - ie z + v k ) 2 

g k (uj)=e k —— — -—j ; 2a; ) \ (3.2) 

T(2v k ) r(i - ie z - v k ) 

where 



1 to 2 + k 2 

Vk = ^l + { Z + l)(z + 2) - e ^ (3 - 3) 
and e z = -, . The result is similar to the RN-AdS case, where the IR Green's function scales with 

z {z+1)k ' 

frequency for arbitrary spatial momenta. 

Mixing between the scalar and the gauge field leads to a negative effective mass squared for the scalar 



4 



> — 



alternative quantization allowed 



H 1 h 



m| F I (z+i)(z+2) m BF +z 2 

4 



Oscillatory region 



Figure 2: Alternative quantization region versus oscillatory region for Lifshitz4 and general dynamical expo- 
nent. These regions never coincide (as the case of RN-AdSs), thus for the case of Lifshitz4 we can always 
start from both alternative and standard quantization fixed points. As z increases, the region of alternative 
quantization grows faster than the oscillatory region. In our case of interest (z = 2), alternative quantization 
is allowed for —4 < m 2 < and the oscillatory region lies in —4 < m 2 < —3 window. 



field. So in a region of the parameter space, fk can become imaginary for large enough values of e z that 
results in possible scalar pair productions and thus instability of the IR geometry. This process is similar 
to what happens in the RN-AdS black-brane case [HHHQIJ]- In contrast with this case in which the gauge 
field couples to matter fields in the near horizon geometry, there exists some examples where the gauge field 
vanishes in the near horizon limit, so there is no pair production and thus the geometry remains stable (an 
example is studied in [2TJ1I2"!"] ) . 

The z dependence of shows that the net effect of increasing the dynamical exponent pushes the system 
toward more stable regions (see Fig. [1]) Q. Also notice that for large values of z, for arbitrary mass and 
charge, the system becomes stable. 

3.2 UV region and matching 

In order to match the solutions of the inner and outer regions at leading order in w-expansion, we have to 
solve the Klein-Gordon equation (|4.12p in the outer region at ui = which readi§ 



r-^dr (r*+ 3 /^ (0 )) + 



(z + 2) 2 e 2 / 1\ 2 2 k 2 

Z 2 r 2 Zf y 1 -^) ~ m -72 



(0) 



0. (3.4) 



In what follows we will restrict ourselves on the z = 2 case where the above equation could be solved exactly. 
Near the boundary (r — > 00) the solution behaves as 



<f>{r ^00)= Ar A - +Br A +, A± = -2 =F vu = -2 T + m 2 , (3.5) 

where the coefficients A and B represent the source and response functions respectively. The scaling dimension 
of the dual operator is real- valued when m 2 > —4, where the vacuum of the dual theory could be stable [7]. 
This bound is called the generalized z-dependent Breitenlohner-Freedman bound which in a d+l-dimensional 
Lifshitz space-time takes the following general form [22] 

2^2 2 (d + z- l) 2 
m > m BF ; m BF = — 2 . (3.6) 

In the case of AdS^+i this bound is m 2 > —j^z, thus the bound is less restrictive on Lifshitz geometries 
compared to AdS ones. In both types of geometries, the bound on the whole space-time does not guarantee 



4 One can see from 112.411 that for z < 1, the gauge field becomes complex, so we always consider z > 1. 
5 From now on we will drop the (0) subscript which indicates the aj-expansion leading order. 
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the stability of the IR region which contains an A0IS2 factor. In our case the BF bound on the near horizon 
geometrjjfl is m 2 > —3 while the whole geometry is stable for m 2 > —4. Thus there exists a window 
(—4 < m 2 < —3) that the four dimensional space is stable while its near horizon geometry becomes unstable. 
The alternative quantization, referred to as Neuman boundary condition is also allowed in the range [23128] 



m| F < m 2 < ™| F + z 2 . 



(3.7) 



In our case this window is —4 < m 2 < which does not coincide with the window of oscillatory region Fig. 

m 



The exact solution of Eq. (13.4[) is 

„2 , 2 



= (r - I! 



+c 2 (- 



■ , - _ „ , 1 + A+ ie 1 + A. 



-)~F(u k 



1 



ie 



ie 2ie r 2 + 2 
~ ¥' ~ ~3~' 3 



ie 2ie r 2 + 2 

T'^IT'^ - 



(3.8) 



where 



= 7^ 



(> 



t 2 + 3fc 2 + 9-4e 2 . 



(3.9) 



Expanding this solution in the near horizon limit and by using Eq. (I4.16[) . one can read c\ and 02.. Considering 
the near boundary limit of the solution p.8[) , one can read the retarded Green's function of the full geometry 
at leading order in uj. The retarded Green's function takes the form of the well-known master formula of [IB] 
which is 



&V 0) (fc) + fo w (fc)s fe M 



(0), 



G R (u,k) 



aV 0) (fc)+a w (fc)S fc (u;) 



(0), 



(3.10) 



where in our case we find 



if aL° )N 
Jo) .(o) 



r(- 



25+"*3 2 ' /fc + 



\r(- 



1+A_ 



IT3— 



-22 ~ v k3~ 



1 + A 4 



-Tr(-2 t / fc )r(-i-A 4 



l+A 4 



+^+f) r(-^-t-^-f )r( ^t- Vh +^) 

i-A_) 2i~" k 3~ 2 " k+: ^2 L r(-2i/ fc )r(-i-A_) 
+"*+#) r(-ii^-y fc -^)r(-ii^-y fc +^") / 



(3.11) 



One can easily check that the determinant of the above matrix is equal to Vk/^Ut which is a consequence of 
equality of the Wronskians in the infinity and the horizon of p.4[) . As long as vu is real- valued, the boundary 
condition of the scalar field in the UV region and thus these coefficients become real. By expanding Eq. 
(|3.10p near u = one finds 



G r {uj -> 0,k) 



af(k) 



(3.12) 



In the case of real IR scaling dimension, the spectral function reads 

v k YvaQk{u) 



A(uj, k) 



v u a { ° )2 (k) 



(3.13) 



^Remember that our AdS2 radius is /„ 



12. 
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thus the near horizon fluctuations are dominant in the low energy regime and the static susceptibility vanishes. 
In the case of complex- valued IR scaling dimension, the matching coefficients become complex and thus the 
static susceptibility bifurcates into the complex plane [12) . 



4 Instabilities, phase transitions and QCPs 

In this section we investigate the analytic behavior of the low energy Green's function which can be used to 
find zero modes (locus of Fermi surfaces) and explore the existence of quasi particles as excitations around 
it. Also one can extract information about the possible quantum phase transitions and its QCPs. In what 
follows we will study the QCPs mentioned in Sec. 1 and try to explain their essences in our model. 

4.1 Zero modes and phase diagrams 

In order to find zero modes of a bosonic system, we must find the poles of the retarded Green's function at 
zero frequency. Following the terminology of |13j . we will refer to the location of these superfluid surfaces 
by ks similar to the location of Fermi surfaces, fop, in fermionic systems. A simple analysis shows that the 
poles are due to the divergences of &!? (&). After some gymnastics with the Gamma functions one finds that 
the superfluid surfaces are located &{j 



where n is a non- negative integer. Note that unlike the case of extremal RN-AdSs black-brane [T3], the 
position of superfluid surfaces does not depend on the charge of scalar field explicitly, although it depends 
implicitly through the IR scaling dimension. The necessary condition for ks to be real requires v>k > which 
implies 



This condition shows that all the zero modes are above the generalized BF bound (|3.6p . In the n — case 



1 + A 
2 



+ u k = -n, 



m 2 > (2n+l) 2 -4. 



(4.1) 



4 = 3(1 + A_) 2 + ^- 



— m 



(4.2) 



where ks vanishes for e = and m 2 = —3. 

The critical value for momentum where the IR scaling dimension becomes imaginary is 



2 e 2 

k c = — ~ m - 3, 



(4.3) 



where 



e 2 >^(m 2 + 3). 



(4.4) 



7 We have used the following identity 1241 




where a and (3 are real functions. 
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Figure 3: Phase diagrams of a charged scalar in Lifshitz background (z — 2) with standard quantization for 
different mass parameters. The solid lines denote the zero modes while the blue region corresponds to the 
oscillatory region. 

Note that we always have ks > k c and also remember that k < k c denotes the oscillatory region [18] . Outside 
the oscillatory region, zero modes must satisfy both conditions described in Eqs. (|4.1I) and (|4.4[) which leads 
to 

e 2 >3(n+l) 2 -J. (4.5) 

These constraints enlightens some features of the parameter space which is shown in Fig. [3] Note that 
according to Eq. (I3.5|) . we always consider m 2 > —4, to ensure that the vacuum of the dual field theory is 
stable. If we deform the system slightly near its zero modes some interesting behavior could occur. If the 
system supports poles in both the upper and lower half of the complex w-plane, it is a sign of instability. In 
order to further study the stability conditions of the system we have focused on the poles in the complex w- 
plane. To do so we have solved the scalar field retarded Green's function numerically to all orders in frequency. 
This could be done by imposing the near boundary solution of the IR region as boundary conditions of the 
UV region. The system supports several poles for different set of parameters. For a given mass parameter, 
at e = and k — 0, the poles lie in the lower half plan while increasing e pushes them upstairs. For a given 
mass parameter and a large enough e, the schematic result is shown in Fig. |4j As |fc| increases, the poles 
approach the normal mode (origin) and by further increasing they finally jump into the lower half plan. 

4.2 Bifurcating quantum critical point 

When the BF bound of the near horizon geometry is violated, the Green's function bifurcates into the 
complex w-plane, known as bifurcating quantum phase transition. At this critical point the effective AdS2 
mass vanishes at zero spatial momentun|f|. At such a critical point the zero frequency Green's function does 
not diverge but bifurcates into the complex plane. The effective AdS2 mass could be tuned by the value of 
the scalar charge. Defining 

v = — y/ 3m 2 + 9 — 4e 2 = y/u, 
6 

a bifurcating QCP happens at 

u = u c = 0. 

8 We set k = 0, in order to neglect the effect of the extra R 2 factor. 
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Figure 4: Schematic plot of the poles of the Green's function for large scalar charge e. This bahavior is found 
by solving the scalar field equation of motion numerically in the UV region, using the near boundary solution 
of the IR region as boundary condition. These poles correspond to the condensed region (u < 0) and the 
most right pole refers to u = 0. 



Near v = 0, the matching matrix could be expanded 



'a™ a'°>\ fa t» N 

5, s r»» I e-"»-l fl „ 



a a 



where 



a 



V2r(-ii^-f)r(-i: 



2 > :S 



a = — - 



P = 



V2T 



1 + A_ 



f)r(- 



1+A_ 



3-^r(^) 




+ i> 




+ 27 -In 18 


V2r( 1+ 2 A + 


f)r(-^ + f 


) 





if) 


{- 


1+A_ ie^ 

2 3 y 








+ 2 7 - In 18 


V2T 


(- 1+ 2 A - 




r 







and ip(x) denotes the digamma function. In this case the determinant condition implies that a/3 — (3a = -^j- 
Expanding the IR Green's function near the bifurcating critical point leads to 



Gk=o{u) = -l + 2i/Q (u>) 

Qoiu) = -H-2iuj). 

Thus the Green's function for the whole geometry near such a critical point takes the form of 

0So(w)+j8 



(4.6) 



G 



aQo(uj) + a 
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Also one can find similar results for the finite temperature case by using Qj. from Eq. (|4.18l) . In this case 
one finds 

-,t i \ t \ , ( 1 is iui \ ( —iui\ 



S » - g o{u) -^[- + J -— + m — . (4.7) 



In the neutral case it reduces to 



gZ(u>) = -7 - fl> Q - ^) - In (ttT) . (4.8) 

As it could be seen in Fig. |4l in the condensed side (u < 0), the system supports infinite number of poles in 
the upper half u>- plane. This is similar to the case of RN-AdS black-brane [T^] . 

4.3 Hybridiezed quantum critical point 

One way to achieve the zero modes (at k = 0) is to add a double trace deformation to the CFT [25^1 . Under 
double trace deformations, the dual theory flows into a new CFT where the scaling dimension of O changes 
from A to d — A. As a result of adding a double trace operator 



2 



d d x O 2 , 



to the action, the Green's function in the momentum space becomes a geometric sum over connected diagrams. 
In this case the Green's function deforms as 

Gf{u,k) ' 



G- 1 (uj,k)+K+' 



where Gr(uj, k) is the two point function for O in the absence of the double trace deformation [3]. 

The double trace deformation displaces the poles of the Green's function and results in new zero modes. 
At the leading order of frequency one finds 



gM {u k) = bV>{k) + b*\k)g k { u ) 

R ' af{k) + bf{k)n+ + (^\k) + b^\k)K + )g k {u})' 
The poles of the deformed Green's function lie at a? (k) + bl (k)K+ — 0, which leads to 



(4.9) 



3 



-vxj 



r(-i - A + )r(-^ + v k - f )r(-^ + Vk 
r(-i - A_)r(-i^± + Vk - |)r(-i^± + Vk 



Hybridized critical points correspond to k+ = k c . Near this critical point, the Green's function at zero 
frequency and momentum limit takes the form 

G [ k\lo -> 0, k -> 0) = ' 



k 2 h k + uhu, + lj 2 Ii w 2 + hC{v)(-~iu) 2v 



3 It could also be achieved by considering superfluid velocity. 
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Figure 5: Phase diagrams for standard and alternative quantizations of a neutral scalar. Left plot (standard 
quantization): For k+ > 0, if there is no n c > 0, the system is unstable (light blue regions). When k c > 
exists, the system is unstable for < k+ < k c (dark blue regions). For n+ < 0, the system becomes unstable 
if k c < exists and k + < k c (red regions). For u < the system is unstable because of the oscillatory region. 
The system is stable in the white regions. Middle plot: is the left plot magnified near the origin. Right plot: 
is the same as the left plot for alternative quantization. 



where 



h k 2 = 



,<<>) 



fc=0 



a (0) 

h = - L - 

" ~~ h (o) 



h u ,2 = 



k=0 

K+ —K c 



~(2) 



h = 



k=0 



~(0) 

a_ 

75T 



fc=0 



vub 



(0) 



2 ! 



and Q = C{v){—iuj) 2v . Since Eq. (|3.4I) is parity invariant, there is no linear term of momentum in the above 
expansion and a similar reasoning leads to the absence of linear frequency term in the neutral case. 

In the case of v < 1, for a neutral scalar the h U 2 term can be neglected and the pole of the Green's 
function originates from the IR region 



hC(u) 



In this case hC(y) > 0. So for k + k c we have Im(w*) ^ 0, which shows that unstable region corresponds 
to k + < K c . 

To further investigate the hybridized critical point, remember that the alternative quantization is also 
allowed in the —4 < m? < window (see Fig. [5]). In Fig. [S]the critical value of k± is plotted as a function 
of u. In the left plot, one can find different critical points in the alternative quantization. The u = line 
refers to the bifurcating critical point while the border specified by k c shows the hybridized critical point. 
The point of intersection corresponds to the marginal critical point. As it was shown, while k+ crosses down 
the k+ — n c curve, the pole of the Green's function crosses the origin of the complex w-plane causing an 
instability. This instability is guaranteed because k c is a single valued function. Obviously this description 
is valid while our analytic calculation is valid in the low energy limit. The middle plot is the same as the 
left one in a wider u region where the zeroes of k c could be seen. The right plot is the same as the left one 
for the case of alternative quantization. Note that in all plots of Fig. [S] we have assumed the UV and IR 
instabilities are related. If they are completely unrelated, there is no stable region for K+ > (see [26] 



3 Wc thank Jie Ren for bringing our attention to this point. 
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4.4 Marginal quantum critical point 

Approaching the bifurcating and hybridized quantum critical points simultaneously, a new critical behavior 
occurs known as marginal quantum critical point. More precisely one must consider both u — > and k + — > k c 
limits at the same time. Since the IR scaling dimension is 5k — \ + i^k, in the Vk — >• limit the double trace 
operator becomes marginal. In the following we will read the retarded Green's function in this limit. 

Consider a neutral scalar on the border of the oscillatory region (to 2 = —3), the solution in the inner 
region can be written in terms of Hankel functions. After considering ingoing boundary condition on the 
horizon one finds 

MS) ~ v^tffTW 

In terms of the UV region coordinates the above solution reads 

MZ -> 0) ~ I 12 ('' - [So + In 12 ( r ~ *)] > 
where Go = —7 — ln(^ £ ). In the outer region the solution is 

1 / r 2 -l\ 

By matching the IR and UV regions one can read the retarded Green's function in this case as follows 

G = ^-21n6 + 6 + 7 -HL . (41Q) 

(k+ + 2)(lnw - 21n6 + 7 - f ) + 6k+ 

The critical point occurs at k + = — 2 (see Fig. [5]). The logarithmic term in Eq. (|4.10[) indicates the bahavior 
of marginal Fermi liquids which is discussed in |12j . 



Conclusions 

In this paper we have studied holographic quantum critical points of a charged Lifshitz space-time. We have 
studied the near horizon geometry showing that it contains an AdS2 factor (as in the RN black-brane case) 
which is dual to conformal quantum mechanics. We considered a charged scalar field in this background and 
worked out the z-dependence of the scaling dimension of the dual operator. In the case of z = 2, we worked 
out the exact solution in the UV and IR regions and computed the analytic retarded Green's function in the 
leading order of frequency expansion. We have found the corresponding zero modes analytically and studied 
the stability regions of the model in the parameter space. We have also studied the bifurcating, hybridized 
and marginal quantum critical points of the system. 



We can summarize the properties of the model as follows: 

• There are several exact information about different quantities in the model such as zero modes and the 
locus of superfluid surfaces. 

• This non-relativistic theory flows to an IR fixed point, but the dynamical exponent still affects the 
bosonic instabilities. The z dependence of the IR scaling dimension of a scalar operator shows that 
strongly correlated systems described by larger z are more stable under scalar condensation. 
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• Stable IR geometry of this model supports both standard and alternative quantizations, unlike the case 
of RN-AdS 5 . 

• The system supports bifurcating phase transition where the generalized BF bound is violated in the IR 
geometry. It also supports a hybridized QCP via adding double trace deformations to the model. 

• By studying these phase transitions one can not distinguish between relativistic and non-relativistic 
field theories in the UV region. 

To further study this background one can probe it by fermionic fields and investigate the existence of Fermi 
surfaces and its QCPs. One can also study scalar instabilities in more general non-relativistic backgrounds, 
such as those with hyperscaling violation. 

Note added: After submitting this paper, T. Andrade and S. Ross brought our attention to [25], where 
although the normalizability seems to allow masses in the range (|3.7[) for alternative quantization, a more 
detailed study reveals that there are severe instabilities unless m 2 BF < m 2 < m BF + 1 . This is essentially a UV 
effect in the field theory and this instability will affect any solution which asymptotically approaches Lifshitz 
solution with the Neumann boundary condition. This in particular implies that in Fig. [5J the alternative 
quantization and oscillatory regions coincide. 
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Appendix 

In this appendix we compute the Green's function in the near horizon region. We do the calculation for a 
scalar and fermion field. These Green's functions in the RN-AdS black-brane case has been done in [27] . 

Scalar field 

Zero temperature case 

We consider the action of a charged scalar minimally coupled to the gauge field as 

S = - J d A x^ [g MN (d M + iqA M ) 0* {d N ~ iqA N ) + m 2 0*0] (4.11) 

so the equation of motion after Fourier transforming like <p(r, x^) = e~ iut+ w <f)(r), becomes 

□0 - {^= d M (V=gA M ) + 2iqA M d M + q 2 A 2 + m 2 ^j = 0, (4.12) 

where the scalar d'Alembertian is defined as □ = ^=9jvf (\/—gd M ) . 

Now consider a scalar field in the near horizon geometry of the Lifshitz space-time Eq. (|2.7j) . Also notice 
that we only consider A^ as the gauge field that can be coupled to the scalar. The equation of motion takes 
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the form of 



e<t>"+ K + e, z f - (m 2 + k 2 )l 2 z 4> = 0, (4.13) 

where e z = -^j-j = ( Z +\) K ■ Near the boundary of AdS2 (i.e. £ — > 0), the solution of the above equation 
becomes 

tftt -> 0) = + , (4.14) 

where cf>o and 0i are the source and the response functions respectively and 



v k = y i + (ra 2 + k 2 )ll - e 2 . (4.15) 

The conformal dimension of the dual operator is 8 k — \ + Vk- The near boundary solution Eq. (|4.14[) . near 
the matching region (£ —> 0), up to an overall factor takes the form of 

4>(£ -> 0) ~ [12(r - + [i2( r - l)]-*-*0 fc . (4.16) 

Eq. (|4.13p could be solved exactly in terms of Whittaker functions as 

0(0 = c in W i&x>Vh {-2iu£) + c out W- ies , Uk (2iu;£), (4.17) 

where applying the ingoing boundary condition at the horizon leads to c ou t — 0. Expanding this solution 
near the boundary of the AdS2 region, one can read the IR Green's function as 

T(—2v k ) T(i - ie z + v k ) „ 
1 (2v k ) T(i - ie z - Vk) 

As in the relativistic case, the IR Green's function scales with frequency for arbitrary momenta. 
Finite temperature case 

In this case we consider the following metric as the near horizon limit at finite tcmprature of the Lifshitz 
black-brane 

ds 2 = gl dT 2 , de , h _ iLr 2 



e \ v a J fi-g N 1 - 



where the gauge field is 

v ~ (z + 1)< ^ e 



A^ = —±—fl « 
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and its Hawking temperaturj^l reads T — (27r^) 1 . The equation of motion for a scalar field in this 
background is 



+ 



24 



1 1 



i « 2 



which its solution is 



Cl i 



F 



v k - ie z , ^ - v k + ie z - iu£ h , 1 - 2v k , 



24 



24 



l 2 4 + 4/i 



and ^fc was defined in (14. 15[) . The ingoing boundary condition leads to 

c _ 2 „ r(-2i/ fc ) T{\ - ie z + Vk)^{\ + ie z + v k - ^60 g 

r(2^ fc ) r(| - ie 2 - v k )Y{\ + ie z - v k - iui£, h ) 

and the retarded Green's function takes the form of 



Y{2v k ) T(\-ie z -v k )T(\ + ie z -v k 



2-kT t 



(4.18) 
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